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A SLIP  MODEL  FOR  FINITE-ELEMENT  PLASTICITY1 

By 

2 

Hendrik  van  Rij 
and 

3 

Philip  G.  Hodge,  Jr. 

Abstract 

Conventional  finite-element  models  are  based  on  displacement 
or  velocity  fields  which  are  at  least  continuous.  However,  it  is 
known  that  perfectly  plastic  materials  may  exhibit  discontinuities 
of  the  tangential  velocity  component  along  certain  lines.  In  this 
investigation,  a two-dimensional  finite  element  model  is 
proposed  which  will  allow  for  such  discontinuities. 

A regular  pattern  of  triangular  elements  is  assumed.  The 
elements  are  assumed  to  be  rigid,  and  across  the  line  separating 
any  two  adjoining  elements  the  normal  displacement  component  is 
continuous,  but  a discontinuity  may  exist  in  the  tangential 
component.  The  defining  equations  - compatibility,  equilibrium, 
and  constitutive  - are  developed  with  the  aid  of  the  Principle 
of  Virtual  Work. 

Prandtl's  punch  problem  is  solved  under  plane  strain  conditions. 
Comparison  is  made  with  existing  analytical  and  other  numerical 
solutions,  in  order  to  evaluate  the  merits  of  allowing  for  dis- 
continuities 

1This  research  was  sponsored  by  the  Office  of  Naval  Research;  the 
results  are  taken  from  a dissertation  submitted  to  the  University 
of  Minnesota  by  one  of  the  authors  (H.v.R.)  in  partial  fulfullment 
of  the  requirements  for  a Doctor  of  Philosophy  degree. 

2 

Engineer,  Sichting  Energieonderzoek  Centrum  Nederland,  The  Netherlands. 
Formerly  Research  Assistant,  University  of  Minnesota 

^Professor  of  Mechanics,  University  of  Minnesota 


1 


1.  Introduction . The  basic  concept  of  the  finite 
element  method  is  that  a given  structure  or  continuum 
is  modeled  by  a finite  assemblage  of  individual 
components  or  elements.  Within  each  element  and  in 
the  reactions  between  the  elements  simplifying 
assumptions  are  made  in  the  static,  kinematic,  or 
constitutive  equations. 

In  a pioneering  paper  in  1956,  Turner,  Clough, 

Martin, and  Topp  [1]  proposed  a finite  element  model 
for  two  dimensional  problems  which  we  shall  refer 
to  as  the  classical  model.  The  elements  are  triangles 
and  the  model  has  the  following  properties:  the 
displacement  field  is  linear  in  each  triangle  and 
continuous  over  the  whole  domain;  the  strain 
stress  fields  are  constant  in  each  triangle;  if  the 
given  body  and  surface  force  field  is  approximated 
by  concentrated  forces  acting  only  at  the  nodes, 
then  these  forces  must  be  in  equilibrium  at  each  node. 

The  method  has  been  extended  to  include  plastic 
behavior . 

Since  the  classical  model  and  the  slip  model 
proposed  herein  are  each  based  on  a field 
which  is  kinematically  admissible  for  the  original 
continuum,  the  results  will  satisfy  any  kinematic 
minimum  principle  for  the  material  involved.  In 
particular,  for  an  elastic/perf ectly-plastic  material 
the  model  yield-point  will  be  an  upper  bound  on  that 

‘Numbers  in  square  brackets  refer  to  items  listed 
in  the  bibliography. 
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of  the  continuum  [2] . However,  as  first  pointed  out 
by  Nagtegaal,  Parks,  and  Rice  [3],  this  upper  bound 
may  be  very  poor  or  even  infinite  if  the  model  does 
not  allow  for  incompressible  motions. 

This  constraint ,e ^ = 0,is  not  easily  satisfied. 
Indeed,  in  a large  field  of  4N  triangles,  there  will 
be  approximately  2N  nodes  and  two  degrees  of  freedom 
per  node  so  that  there  is  an  average  of  one  degree 
of  freedom  per  triangle.  Therefore,  one  incompres- 
sibility constraint  per  triangle  will  leave  essentially 
no  degrees  of  freedom  in  the  total  domain,  regardless 
of  the  number  of  elements. 

Nagtegaal  et  al  resolved  this  dilemma  by  re- 
quiring the  triangles  to  be  grouped  in  quadrilaterals 
such  as  the  squares  in  Fig.  1.  It  is  easily  shown 
that  if  three  of  the  four  triangles  in  a quadrilateral 
undergo  any  incompressible  deformation,  the  fourth 
will  necessarily  satisfy  e ^ = 0.  Thus  there  are 
only  3N  constraints,  which  leaves  one  degree  of 
freedom  for  each  four  triangles. 

The  velocity  field  associated  with  the  yield- 
point  load  of  many  continuum  problems  exhibits 
another  phenomenon  which  is  not  well  modelled  by  the 
classical  finite  element  model,  namely,  discontinuity 
of  certain  tangential  velocity  components.  For 
example,  let  us  consider  the  identation  of  a semi- 
infinite domain  by  a rigid  beam  under  conditions  of 
plane  strain.  For  the  simple  case  of  a rigid/perfectly 
plastic  material,  Prandtl  [4,2]  presented  a yield- 
point  incipient  flow  field  (see  Fig.  8b)  which 


3 


‘ 

I 


involved  part  of  the  domain  (the  triangles)  moving  as 
rigid  bodies,  and  part  of  it  (the  quarter  circle) 
deforming  with  zero  volume  change.  We  observe  that 
the  deforming  part  of  the  material  moves  relative  to 
the  non-deforming  part  with  a finite  tangential 
velocity  discontinuity  across  the  boundary  separating 
the  regions. 

Therefore,  in  view  of  the  facts  that  any  continuous 
finite  element  model  has  difficulty  representing 
incompressible  flow,  and  that  tangential  velocity 
discontinuities  actually  exist  in  continuum  solutions, 
it  appears  worthwhile  to  investigate  a finite  element 
model  which  will  allow  for  such  discontinuities. 

• In  this  paper  we  present  a slip  model  in  which 
these  discontinuities  are  the  only  permissible  motion, 
i.e.  we  consider  an  array  of  triangular  elements  which 
are  internally  rigid,  but  can  slip  relative  to  each 
other.  In  Secs.  2 and  3 we  derive  the  basic  equations 
in  the  interior  and  on  the  boundary  based  on  the 
Principle  of  Virtual  Work.  In  doing  this,  we  visualize 
the  edge  between  two  triangles  as  a rectangle  of  small 
thickness  and  ignore  second-order  effects.  Also,  it 
is  shown  that  we  have  one  degree  of  freedom  per  node 
for  this  incompressible  model.  Section  4 considers  a 
simple  problem  explicitly  to  point  out  an  interesting 
aspect  of  non-uniqueness,  and  in  Sec.  5 we  present 
results  for  an  approximation  to  the  Prandtl  problem 
for  a rough  punch.  The  results  are  compared  with 
known  exact  and  numerical  solutions.  Finally,  some 
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conclusions  are  drawn  in  Sec.  6, 


2.  Interior  domain.  Although  the  theory  can  be 
developed  for  any  regular  or  irregular  arrangement 
of  triangles,  we  shall  consider  only  the  case  of  a 
regular  array  of  isosceles  right  triangles  illustrated 
in  Fig.  1.  In  this  two-dimensional  finite  element  model 
the  triangles  are  considered  as  rigid  bodies  which  can 
slip  relative  to  each  other  but  must  not  become 
separated.  Therefore,  the  only  "strains"  in  the  model 
will  be  measured  by  the  relative  motion  or  slip 
between  the  two  adjacent  triangles. 

Let  d^j  be  the  slip  along  the  edge  ij  joining 

nodes  i and  j.  It  will  be  taken  as  positive  when, 

from  the  viewpoint  of  an  observer  at  node  i and  facing 

edge  i j , the  triangle  on  the  right  moves  towards  the 

observer  with  respect  to  the  triangle  on  the  left. 

Evidently  d..  = d...  A dimensionless  strain  is  then 
ID  31 

defined  by 

a).  . = d.  . 1.  ./I  2 (1) 

13  ID  ID  0 

where  is  the  length  of  the  edge  and  lQ  is  a 

reference  length  for  the  domain. 

If  we  ignore  boundary  effects,  then  in  a large 
domain  with  N nodes,  there  will  be  2N  triangles  and 
3N  edges,  hence  there  will  be  3N  generalized  strains. 

The  requirement  that  triangles  do  not  separate  will 
introduce  several  compatibility  constraints  among 
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these  strains.  In  the  first  place,  there  can  ob- 
viously be  no  relative  rotation  of  triangles  so  that 
at  most  we  need  consider  only  two  degrees  of  trans- 
lational freedom  per  triangle,  or  a total  of  4N  for 
the  field.  However,  since  even  a uniform  separation 
is  prohibited,  there  will  be  an  additional  constraint 
on  each  of  the  3N  edges,  so  that,  we  are  left  with  only 
N degrees  of  freedom,  i.e.  one  degree  of  freedom  per 
node.  It  follows  that,  if  one  elementary  mechanism 
is  associated  with  each  node,  all  possible  infinitesimal 
motions  in  a large  field  may  be  expressed  as  linear 
combinations  of  these  mechanisms. 

Figure  2 shows  the  elementary  mechanisms 
associated  with  generic  "small"  and  "large"  nodes. 
Evidently  Fig.  2b  may  be  obtained  from  Fig.  2a  by 
a change  of  scale  and  a rotation  so  that  the  mechanisms 
are  essentially  the  same. 

We  define  generalized  displacements  by 


eD  = id/i0‘ 


e = ia/i' 

G o 


(2) 


The  displacements  across  AD  and  DG  must  be  equal  so 
that  we  can  refer  directly  to  the  diagonal  edge  AG 
in  Fig.  2a;  the  diagonals  through  G,  however,  must 
be  considered  independently.  Evidently  the  strains 
associated  with  the  mechanisms  in  Figs.  2a  and  2b 
respectively  are 

^FG  = WGB  = WBA  = U,AF  °D 
WAG=WBF  = - 2°D 


(3a) 


~ ’ r '~"~r  r ' r 

' * 


/ 


“fl  ^lh  whb  wbf  gb 


WFG  WLG  WHG  WBG  °B 


(3b) 


The  total  strain  along  any  edge  will  be  the 
resultant  of  several  mechanisms.  For  example 

(4a) 


to  =0  + 0 - 0 — 0 

FG  D I F G 


“AG  - eB  + 9P  - 26D 


(4b) 


Let  us  expand  a generic  edge  PQ  to  a finite 
thickness  6 as  in  Fig.  3.  Then  the  displacement 
field 


u = (1/2) (u’  + u”)  + y dpQ/£Q6  v = V (5) 

is  continuous  with  the  rigid  triangle  motions  and 
the  internal  strain  is  the  constant 


Yxy  dPc/*!'o<5  ^Pc/c/^PQ 5 


(6) 


Since  a constant  strain  implies  constant  stress  for 
any  reasonable  constitutive  behavior,  we  define  a 
generalized  stress  for  the  edge  by 


tPQ  V1™  k £_£ 6 


ct  dA  a 
xy  _ xy 

“ k 


(7) 


<S+0  A o PQ 
where  k is  the  yield  stress  in  shear.  Then  the 
dimensionless  internal  work  is 


w . , = W . . /k£  = Trirt(jo_- 

int  int'  o PQ  PQ 


(8) 


which  satisfies  Prager's  criterion  [ 5]  for  generalized 
variables . 

No  external  work  will  be  done  during  the 


A 


mechanism  motions  in  Fig.  2,  so  that  the  associated 
internal  work  must  vanish.  Thus  we  obtain  the  generic 
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equilibrium  equations 


TFG  + TGB  + TBA  + taf  2tag  2tbf  ° 


(9) 


tfl  + tlh  + thb  + tbf  ” tfg  " tlg  thg  tbg  0 


at  small  and  large  nodes,  respectively. 

Since  stress  and  strain  are  constant  in  an  edge, 
the  constitutive  equations  for  an  elastic/perfectly 
plastic  material  may  be  written  as 


IF  1 T^l  < 

1 OR  T_ _ _ < 0 

1 PQ 1 

PQ  PQ 

THEN  XpQ  = 

<g'V“pq>“pq 

(10a) 

ELSE  ipQ  = 

0 

(10b) 

where  we  have  defined 

G'=G/6  (11) 

Equations  10  and  11  imply  that  the  edge  is  made  of 
a fictitious  material  whose  shear  modulus  tends  to 
zero  with  6 so  that  G'  remains  finite.  For  such  a 
model  non-zero  slips  can  be  found  even  in  the  elastic 
range,  but  they  will  be  expressed  in  terms  of  the 
single  unknowable  modulus  G' . 


3.  Boundary  value  problem.  In  the  previous  section 
we  have  defined  a model  with  one  generalized  displacement 


variable  6^  for  each  of  the  N interior  nodes.  For  N 


sufficiently  large,  there  are  approximately  N generalized 
displacements  0^,  2N  generalized  strains  uk  ^ , and  2N 


stresses  x...  To  determine  these  5N  unknowns , we  have 
il 


2N  strain-displacement  equations  (3) , N equilibrium 
equations  (9) , and  2N  constitutive  equations  (10)  . 
However,  for  any  finite  number  N,  the  number  of 
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variables  and  equations  may  vary  slightly  from  these 
numbers  due  to  some  nodes  and  edges  being  on  the 
boundaries . 

To  illustrate,  let  us  regard  ACHF  in  Fig.  2b  as  a complete 
domain  subject  to  the  indicated  boundary  conditions. 

We  regard  the  boundary  conditions  as  being  applied  to 
fictitious  triangles  external  to  the  boundary.  For 
example,  triangle  E'BC  above  line  BC  moves  down  with 
a prescribed  velocity  as  indicated  in  Fig.  4. 

No  mechanism  is  associated  with  thf  fictitious 
point  E',  and  the  mechanism  associated  ■'ith  B (Fig.  4) 
is  simpler  than  that  for  an  interior  node.  Therefore, 

Eq.  (4a)  should  be  replaced  by 

“BC  =V<eB  + ec)/2  <12) 

Next,  we  observe  that  in  any  finite  rectangle 
the  total  number  of  degrees  of  freedom  is  one  less 
than  the  total  number  of  nodes.  Indeed,  if  each 
interior  and  boundary  node  is  given  an  equal  motion 
0,  the  result  in  view  of  (4)  and  (12)  will  be  to 
leave  the  entire  domain  in  its  original  position. 

Therefore,  we  may  arbitrarily  set  9^  = 0 and  eliminate 
the  corresponding  equilibrium  equation  as  being  a 
linear  combination  of  the  other  equilibrium  equations. 

In  the  interior  of  Fig.  2b  there  are  5 unknown 
0.  and  12  unknown  w . . and  t. . with  a corresponding  5 
equilibrium  and  12  constitutive  and  strain-displacement 
relations.  We  now  examine  each  boundary  segment  to 
determine  the  additional  unknowns  and  equations . 


v- 


- - --  . . '■*  - -d* 


The  fictitious  triangles  J'ML  and  I'LK  cannot 
move  vertically,  so  that  0M  = 0 implies  that  0=0 
which  in  turn  implies  0^=0.  The  stress  xLM  = xRL  = 0 
and  the  strains  are  free  but  of  no  interest.  Thus  we 
have  added  neither  unknowns  nor  equations  to  the  list 
of  unknowns.  Evidently  the  same  conclusion  applies 
to  CHM  and  0 = 0 = 0. 

n C 

On  KFA,  a non-zero  normal  displacement  is  prescribed. 
Since  0 = 0 it  follows  from  Fig.  4 that 

x\ 


0„  = 2u^  i/l 
F o o 


0.  = 2ew  = 4u  £/* 

A F o o 


A mechanism  motion  compatible  with  (13)  leads  to  the 
equilibrium  equation 

TKG  + TAG  " TFG  + 2tFB  = T2  (14) 

where  T^  is  the  averaye  dimensionless  normal  force 
on  AK.  Since  x = x__  = 0 and  the  strains  are 

Kr  x1  A 

irrelevant,  we  have  added  one  unknown  force  T£  and 

one  new  equation  to  our  system. 

On  the  top  edge  x,D  = 0 and  cu  n is  irrelevant, 

AB  Ao 

but  xBC  and  toBC  are  additional  unknowns  subject  to 
(10)  and  (12)  (with  0C  = 0) . The  mechanism  at  B must 


satisfy 


eB  - - 2vol/lo 


and  the  corresponding  equilibrium  equation  is 

2(lAG  + TGC  “ TBG)  “ TBC  = T1  (16) 

where  T^  is  the  average  dimensionless  normal  force 
on  BC.  Therefore,  the  total  system  balances  with 
33  equations  and  unknowns. 
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4.  Solution.  As  with  most  elastic-plastic  problems, 

the  general  method  of  solution  is  to  assume  that  all 

0.,  t.  .,and  a).  . are  known  at  a generic  time  t , 

convert  all  equations  to  rate  form, and  solve  for  0^, 

x..,and  a)..,  and  then  find  0.,  x . . , and  u.  . at  t + At 
1}  13  l 13  13  o 

by  linear  integration.  The  interval  At  is  chosen  so 
that  no  edge  changes  between  elastic  and  plastic 
during  the  interval,  hence  each  edge  is  always  on 
the  same  branch  of  (10)  during  each  rate  solution. 

Problems  associated  with  choice  of  branch  in 
(10)  turn  out  to  be  easily  handled,  but  an  interesting 
lac!;  of  uniqueness  occurs  at  certain  stages.  To 
illustrate  this  phenomenon  let  us  consider  in  detail 
the  problem  where  vq  = 4d,  uq  = -d,  and  we  are  not 
concerned  with  the  resultant  forces  T^  and  Then 

the  non-zero  boundary  rates  are 


0fi  = -4d 


0,  = -2d 
A 


eF  = -d 


(17) 


where  we  have  taken  SL  =21. 

o 

There  are  5 homogeneous  equilibrium  rate  equations 


^FG  + 

^LG 

- 2 

(^KG  + 

^FL} 

= 

0 

(18a) 

^FG  + 

"gb 

- 2 

(tfb  + 

"ag} 

= 

0 

(18b) 

*FL 

+ *LH  + *HB  + 

• 

x 

FB 

_(tfg  + 

^LG 

+ 

thg  + tbg) 

= 0 

(18c) 

’LG  + 

"gh 

- 2 

(*LH  + 

tgm^ 
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0 

( 18d) 

*BG  + *GH 

• 

+ T 

BC 

- 2<;gc 

+ "bh1 

= 0 

(18e) 

and 

13  strain- 

displacement  relations 

“fl  = "2§ 

I + 

»G 

WKG 

= -26 

I 

- a 

(19a, 

unloading  in  some  edges.  Indeed,  if  the  motion  about 
node  J were  the  only  contribution  to  edges  LG  and  GH 
for  example,  0 would  have  to  be  zero.  The  freedom 

• U 

allowed  by  (23)  is  due  to  the  fact  that  other  mechanisms 
also  produce  strain  rates  and  that  only  the  total  strain 
rate  is  subjected  to  the  inequalities  (22)  . 

5.  Example . A computer  program  was  written  to 
implement  the  finite  element  slip  model,  and  was 


used  to  solve  the  problen  illustrated  in  Fig.  1.  An 
elastic/perfectly  plastic  material  is  placed  in  a 
perfectly  lubricated  box  and  indented  with  a rigid 
punch.  This  problen  is  a finite  domain  approximation 
to  the  Prandtl  problem  ft , 2]  of  a rigid  rough  punch 


indenting  a semi-infinite  perf ectly-plastic  material 
under  conditions  of  plane  strain.  It  is  also  essentially 
the  one  considered  in  Sec.  3 except  that  the  domain  is 
larger  and  uq  = 0 . Details  of  the  solution  and  a 
discussion  of  the  computer  program  may  be  found  in  [ 6 1 • 


L3 


The  first  stage  of  our  solution  is  the  fully 
elastic  range.  This  stage  is  terminated  when  edge 
pq  of  Fig.  1 reaches  positive  yield  for  vq  = 2.27  G'/k 
and  an  averaged  compressive  stress  under  the  punch  of 
2.50,  i.e.,  T^  = -2.50.  Figure  5 pictures  the  elastic 
displacement  field  by  showing  the  averaged  vertex 
displacements  at  representative  small  nodes. 

As  v is  increased  above  its  maximum  elastic 
o 

value  there  will  be  a succession  of  elastic-plastic 
stages.  The  sequence  in  which  edges  become  plastic 
is  indicated  in  Fig.  6,  and  the  load-deflection  curve 
for  the  punch  is  shown  in  Fig.  7. 

No  difficulties  with  the  solution  are  encountered 
until  stage  10.  The  first  rate  solution  at  this  stage 
shows  that  edge  pv  which  reached  positive  yield  at  the 
end  of  stage  6,  now  has  a negative  strain  rate  and 

violates  the  inequality  of  the  constitutive  equation 
(10) . Therefore,  this  rate  solution  must  be  rejected 
and  edge  pv  assumed  to  be  elastic  again  for  stage  10. 

It  remains  elastic  during  all  subsequent  stages.  The 
maximum  unloading  occurs  at  the  end  of  stage  15  and  is 
about  2%  less  than  the  yield  stress. 

Edge  hi  alternates  being  elastic  and  plastic.  At 
the  end  of  stage  10  this  edge  reaches  positive  yield 
and  remains  at  yield  until  the  beginning  of  stage  15. 

It  is  elastic  in  stages  15  and  16,  plastic  in  17, 
elastic  again  in  stages  18-21,  and  plastic  from  then  on 
After  initial  yield,  the  elastic  stress  is  never  less 
than  99%  of  the  yield  stress  so  that,  in  fact,  no 
serious  error  would  be  introduced  if  unloading  were 


neglected. 

At  the  end  of  stage  19  we  see  in  Fig.  6 that  the 
edges  of  the  elementary  large  node  mechanism  of  node 
p are  all  at  yield.  Therefore,  0^  is  undetermined, 
subject  to  the  inequalities  of  Eq.  (10) . These 
inequalities  reduce  to 

-. 179vq  > 9p  > - . 189v0  (24) 

for  stages  20  and  21. 

Fig.  6 also  shows  that  at  the  end  of  stage  21 

all  of  the  edges  associated  with  the  elementary 

mechanism  of  small  node  k are  also  at  yield  so  that 

there  is  a further  singularity  in  the  stiffness  matrix, 

The  rate  solution  for  stage  22  puts  the  following 

bouhds  on  6,  and  3 : 

k p 

-.114Vq  0k  >_  --.121vq 


• 055Vq  + 20k  <_  0 <_  -.190vq 

At  the  end  of  stage  22,  four  edges  simultaneously 
reach  yield,  T^  = -6.00,  and  vQ  = 25.05  G'/k.  It 
is  evident  from  Fig.  6 that  sufficient  edges  have 
yielded  to  permit  the  mechanism  motions  shown  in 
Fig.  8a  with  all  inequalities  in  (10)  being  satisfied, 
as  well  as  numerous  other  mechanisms.  Since  all 
stresses  in  the  non-homogeneous  equilibrium  equation 
corresponding  to  (16)  are  at  yield  and  hence  constant, 
the  load  is  constant,  and  T^  = -6.00  represents  the 
yield-point  load. 

In  order  to  evaluate  this  method  of  velocity- 
discontinuities  we  compare  the  results  in  the  fully 


15 


elastic  range  and  at  the  limit  load  condition  with 
analytical  and  other  numerical  results. 


Green  and  Zerna  [ 7 ] have  obtained  an  analytical 
solution  to  the  Prandtl  problem  for  an  incompressible 
elastic  material.  Their  displacements  are  shown  by 
dashed  arrows  in  Fig.  5.  The  dotted  arrows  show  the 
solution  of  the  elastic  problem  according  to  the 
classical  finite  element  model  with  constant-strain 
triangles  and  no  slip,  taking  v = 0.49. 

Figure  5 shows  surprisingly  good  agreement  for 
the  elastic  solution.  It  must,  of  course,  be  emphasized 
that  the  dimensionless  displacements  for  the  other 
solutions  are  defined  in  terms  of  known  physical 
constants,  whereas  the  slip  model  they  include  the 
essentially  unknowable  modulus  G'.  However,  by 
choosing  the  constant  G'  = 14.05  G so  that  the  over- 
all slopes  match  in  Fig.  7,  we  obtain  good  agreement  over 
the  entire  field. 

The  agreement  of  the  three  solutions  is  quite 

good  near  the  punch  and  along  most  of  the  free 
boundary.  The  poor  agreement  along  the  edges  suggests 

that  the  size  of  the  box  in  the  model  was  too  small 

to  give  a good  approximation  to  a semi-infinite  domain. 

Little  information  is  available  for  comparison 
of  the  various  elastic-plastic  stages,  but  a partial 
solution  at  the  yield-point  load  was  first  given  by 
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Prandtl  [ 4 ] and  was  later  shown  by  Shield  and  Drucker 
to  be  complete  [ 8 ] . Fig.  8b  shows  the  resulting 
kinematic  pattern.  The  yield  point  stress  is: 

T1  = - (2  + tt)  = - 5.14  (26) 

compared  with  our  value  of  -6.00. 

6.  Conclusions . Other  examples  have  been  considered 
in  [ 6 ] . These  include  the  Prandtl  problem  for  a 
smooth  punch,  and  tension  of  bars  with  a deep  perpendicular 
slit  in  either  or  both  sides.  The  results  are  all 
generally  similar  to  the  problem  of  Sec.  5.  In  all 
cases,  the  program  ran  smoothly  and  the  difficulties 
posed  by  unloading  and  by  non-unique  velocity  solutions 
proved  easily  surmountable. 

We  have  presented  a model  based  on  slip  which  is 
easily  adapted  to  computer  use  and  which  appears  to 
give  reasonable  results  for  an  elastic/perfectly-plastic 
material.  Potential  difficulties  caused  by  a partial 
lack  of  uniqueness  prior  to  yielding  do  not  prove 
unsurmountable , although  it  would  be  desirable  to  further 
investigate  this  phenomenon.  The  fact  that  G'  cannot 
be  determined  means  that  only  relative  displacements 
can  be  determined,  but  it  does  not  affect  the  stress 
solution  and  appears  a small  price  to  pay  for  the 
simplicity  of  the  model. 

An  obvious  question  arises  as  to  the  feasibility 
and  desirability  of  a model  which  would  include  both 
constant-strain  deformations  of  the  triangles  and  slip 
between  them.  It  turns  out  that  such  a combined  model 
separates  into  two  problems  corresponding  to  its  two 
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constituients , and  that  the  two  problems  are  only 
weakly  coupled.  Details  of  this  combined  model  are 
reported  an  elsewhere  [ 6 ] . 
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